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Title: Seismic Detection Apparatus and Related Method 
Field of Invention 

This invention relates to seismic detection apparatus and to a method of analysing data acquired by 
such apparatus. 

Background to the Invention 

In seismic detection, a seismic source signal propagates tlirough different rock substrates or strata 
within the earth, so producing compression waves (P-waves) and shear waves (S-vvaves) which can 
be analysed to determine the direction and extent of geological features in strata. 

Generally a seismic source signal is produced, either on land, sea, or in a borehole, the source signal 
usually being produced as either acoustic or elastic energy. Some of the energy radiates downward 
through geological layers within the earth and is reflected in varying proportions at different layer 
boundaries and this reflected energy can be detected at the surface. Some of the energy within the 
seismic source signal propagates directly along the ground producing a wave noise signal known as 
ground-roll. To detect the acoustic reflections which are in the form of P-waves and S-waves, a 
linear or area! array of geophones is used to reinforce the reflected energy and to attenuate the 
energy associated with ground-roll so that a good signal to noise ratio can be achieved for the 
reflected wave components of interest. 

It is desirable to be able to separate the P-wave and S-wave components of an elastic wave field as 
the separated components can be used more effectively to identify different characteristics than 
when combined. It is an aim of the present invention to provide improved separation of seismic 
signals into P-wave and S-waves. 

Summary of the invention 

According to one aspect of the present invention, there is provided seismic detection apparatus 
comprising seismic detection means capable of detecting a plurality of seismic components over a 
selected volume. 

Preferably the seismic detection means comprises a plurality of receivers spaced apart to enclose the 
selected volume. The receivers may include hydrophones, geophones, accelerometers or a 
combination of these, and preferably are adapted to detect seismic signals along tliree mutually 
orthogonal axes. Thus the receivers may be provided by three-component geophones. 

Preferably the selected volume is formed by at least tliree receivers spaced within a plane and at 
least one receiver placed external to the plane, thereby to enclose the selected volume over which 
measurement of seismic signals occurs. 

Preferably the receivers are spaced to form a tetrahedron with three receivers spaced to define a 
triangle within the plane and a fourth receiver placed external the plane but in line with the centroid 
of the triangle. Other configurations of spaced apart receivers may be used depending on the 
characteristics of.the wavefield being measured, for example receivers may be spaced to form an 
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octahedral volume, a cubic volume, or a spherical volume. Certain shaped volumes may be more 

appropriate for diFferenTlIsesrfor'e"xahiplFr 

more appropriate for use with permanent downhole sensors. 



The seismic detection apparatus in accordance with the present invention is suitable for use on the 
surface, with the enclosed volume being achieved by placing one receiver further below the surface 
of the other receivers, for example by burying it beneath the surface. With marine and ocean 
bottom surveys, the volume may be defined by tethering the receivers at different water depths, for 
example by placing three receivers spaced apart at the same depth, and a fourth receiver at a 
different depth. The seismic detection apparatus according to the present invention is also of use in 
seismic detection is a transitions zone such as swamps and marshes, where again tethering of the 
receivers at different depths will define the enclosed volume. 

The apparatus may also have all the receivers incorporated into one body for attaching to a w^ireline 
and using downhole, or for using permanently downhole. Such a permanent sensor is of particular 
advantage where monitoring of a reservoir or other geological feature over time is required. This 
may be appropriate when monitoring fluid flow occurring within a reservoir as permanent sensors 
can be used to detect natural seismicity resulting from fluid opening and closing micro-lissures 
within the substrates and producing micro-earthquakes. Whilst these micro-earthquakes produce 
very weak broadband signals, by use of permanent sensors, this activity can be detected and 
production adapted by remedial action such as acid to remove carbonate deposits, or fracturing rock 
with high pressure steam or water, so as to improve production efficiency. 

Preferably the apparatus further comprises processing means which analyses detected seismic 
components from individual receivers to separate P-vvave components from S-vvave components. 
The processing means may be provided at substantially the same location as the receivers to allow 
for on-site processing, or remote from the receivers to allow off-site processing. 

In surveys of the vertical component of seismic data, the receivers may be adapted to detect 
components along one axis and orientation of the receivers selected so as to permit processing of 
one-component data to identify P-wave components. Therefore the receivers may be provided by 
one-component geophones. 

Preferably the spacing of the receivers is selected to be smaller than the wavelength of the detected 
seismic components. A preferred spacing is between 0.05m and 0.50m, although larger spacings of 
up to a third, of the shortest detected wavelength are possible, and thus the range may be between 
1 .5 and 15m. 



In use, the apparatus may be operated to detect seismic signals over a time period of the order of 
five seconds where P-wave components arc of principal interest, with any marginal S-wavc 
component detected at that time being used to produce improvements in the signal to noise ratio. 

Where a detection time of the order 10 to 15 seconds is used, both P-wave and S-wave components 
will be detected within the seismic signals and may be separated for multi-component- inversion 
yielding images of reservoirs and other geological formations, and full reservoir characterisation. 

In accordance with another aspect of the invention, there is provided a method of processing 
seismic data comprising acquiring seismic data relating to a wavefield over a selected volume of 
acquisition, and measuring the curl and divergence of the wavefield from the seismic data, to 
thereby identify seismic components within the seismic data. 



As P-waves are curl free and S-vvaves divergence free, this allows the separate identification of P- 
wave and S-vvave components. Up-going and down-going wavefield components in boreholes may 
also be identified from the seismic data. Attenuation of unwanted seismic components from data 
acquired downliole is also enabled, and noise sources such as ground roll, airwave and inter-sensor 
perturbations can be removed from the data by signal processing. 

The method preferably further comprises averaging the curl and divergence over the selected 
volume of acquisition, and where appropriate acquiring seismic data over a selected period of time. 

Additionally the method may further comprise defining the selected volunie of acquisition to be 
small compared to the wavelength of seismic waves to be detected. 

The invention will now be described with reference to the accompanying drawings in which: 

Figure 1 shows an example of the distribution of seismic data over a specified period of time; 
Figure 2 shows one embodiment of seismic detection apparatus according to the present invention; 
Figure 3 shows a second embodiment of seismic detection apparatus according to the present 
invention; and 

Figure 4 shows a tetrahedron used in explanation of a processing method according to the present 
invention. 

Description 

In Figure 1, a typical hodogram of a three component marine shot gather survey is shown, where 
the X axis is the offset in metres and the y axis is time in seconds from the generation of a source 
signal at t=0. A hodogram displays the different seismic components detected over a specified 
time-w^indow. The first arrivals around t== Is predominantly comprise P-waves 10,10'. Near- 
surface scattering and channelling in slower near-surface layers are revealed in the P-wave arrivals 
contained in the hodogram through some degree of elliptical motion. Following the P-wave arrivals 
is a relatively quiet section 12, 12' (where deeper reflections can be observed) until the S-wave 
arrivals 14, 14' are observed near t=3s. The S-waves also scatter and become channelled in the 
near-surface which causes ellipticity in the hodogram. 

Following the S-wave arrivals are elliptically polarised Rayleigh waves 16, 16' (ground-roll). The 
innermost "cone" 20 of the section consists of strongly scattered energy that typically masks deeper 
reflections. 

Noise in seismic"data io a large degree' comprises source-generated ground-roll and S-waves. 
Separating P-wave data from S-wave data is useful to reduce noise attenuation. Also P-wave data 
may be of primary interest with the complete removal of S-waves being preferred, for example in 
shallow seismic imaging for engineering or environmental purposes. 

Measuring the P-wave and S-wave components of an elastic wavefield separately has several 
advantages. First, separation of P-waves and S-waves enables better and more accurate analysis of 
each separated component thus improving moveout analysis and Amplitude Versus Offset (AVO) 
analysis. Secondly heterogeneity in the near-surface (as well as deeper down in the crust) and 
rough surface topography, scatter S-waves to a greater extent than P-waves, with scattered S-waves 
significantly contributing to noise in seismic reflection data. 
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Figure 2 shows one embodiment of a seismic detection apparatus according to the present 

invention, with 'foITr tlTree-c'oniponent~detectors7Such' as~geophones-or"h^^^ 

form an enclosed tetrahedral vokime. Geophones are velocity motion sensors, with hydrophones^' 
being pressure sensors. Geophones 22, 24, 26 are arranged in a horizontal plane 30 with each 
geophone being spaced the same distance from the other two geophones, so as to form an 
equilateral triangle. The fourth geophone 32 is placed vertically below the centroid 34 of the 
equilateral triangle formed by geophones 22, 24, 26. Placing of geophone 32 external to the 
horizontal plane 30 is achieved by positioning the horizontal plane 30 on or close to the Earth's 
surface with the fourth geophone 32 dug to a position well below surface. Within a marine, ocean 
bottom or swamp environment, the arrangement of receivers is achieved by tethering three 
hydrophones and/or geophones at one depth and one hydrophone and/or geophone at a different 
depth. 

An ideal tetrahedral volume is not essential, and thus for example two geophones placed external to 
horizontal plane 30 can be used. 

The geophones 22, 24, 26 are positioned one metre apart, although spacing will vary depending on 
the surrounding rock substrate and how fast waves travel within the substrate. The spacing of the 
geophones to form an enclosed volume 36 is chosen so that the wavefield is detected on a scale 
which is smaller than the wavelengths of interest. For wave frequencies of between 50 to 100 Hz 
with velocity ranging from 300-2000 m/s, the wavelength tends to be in the range 5 to 50 m. To 
record the wavefield over such a range, the spacing of the geophones is generally a third of the 
wavelength and generally in the range 1.0m to t5m, although may be as small as 0.05 m. The 
tetrahedron arrangement shown in Figure 2 uses the minimum number of geophones required to 
define an enclosed volume, although other configurations with more geophones are possible, such 
as those forming an octahedron, cube or sphere. 

The use of four three-component geophones distributed as shown in Figure 2 allows the separation 
of the wavefield into P-waves and S-waves. After the seismic data has been recorded either where 
acquisition of data is occurring or remote from the geophones, differences in geophone coupling 
can be removed by normalising separately the responses of each component at the corners of die 
tetrahedron. If required, more accurate compensation for coupling differences can be obtained by 
using information from adjacent tetrahedral arrays of geophones, while preserving the offset 
varying character of the seismic data. The geophone coupling may also be frequency dependent 
varying from geophone to geophone. 

The three-component geophones 22, 24, 26, 32 detect components of the seismic signal at tliree 
orthogonal axes. By measuring the curl of a displacement field u (x, t) within the enclosed volume 
36 and the divergence of the displacement field within the volume 36, it is possible to separate out 
the P-wave and S-wave components as P-waves have zero curl and S-waves have zero divergence. 

With the array of geophones shown in Figure 3, the four geophones 40, 42, 44, 46 are one- 
component geophones. The array encloses a tetrahedral volume 50 once more, but the ax is along 
which each one-component geophone measures is calculated and precisely implemented .when 
positioning the apparatus to ensure that the divergence can be deduced from the single components 
measured. Measuring with single component geophones provides less information about the 
seismic signals. However at present 90% of seismic surveys only monitor the vertical component, 
and the apparatus of Figure 3 when used with a listening time of five seconds allows a substantial 
improvement in the signal to noise ratio if looking principally for P-wave components. 
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Once data arc acquired using llic gcophonc army enclosing a volume, separation (^ftiic wnvc-dcld 
into F- waves and S-waves or improvement of the signal to noise ratio Ibr P-wave components, can 
be undertaken by use of the method described below. 

The clastic wave equation for a displacement Field u(.v,/) in an isotropic medium with I.ame 
parameters A and /i and density p can be written as: 

pu = f+(/L + 2//) V(V • u) -// V X (V X u) (1) 

Here f denotes a distribution of body forces. Lame's thereon (Aki and Richards, 1980) states that 
there exist potentials O and T of u with the following properties: 

u = VcD + V X ^F, (2) 

V-4^ = 0, (3) 

O 

CD + a^V^cD, (4) 
P 

^ = ^^J3'V'^V, (5) 

where a and P are the P- and S- velocities. An elastic wavefield u can thus be decomposed into its 
P- and S-wave components, VO and V x 4^, respectively. Equations (2) and (3) yield 

V-u=V-0, (6) 

Vxu = VxVxT (7) 

(Rade, L., and Westergren, N., 1988, BETA p Mathematics Handbook: Student litteratur, Lund.). 
By measuring the curl and the divergence of an elastic wavefield we can thus measure the P- and S- 
wave components separately. 

The technique for P- and S-wave separation is based on recording the elastic wavefield around a 
small volume as compared to the wavelength. The volume may be of arbitrary shape, but the 
equations for wavefield separation are simplified if the volume has a known deterministic shape 
such as that of a tetrahedron, cube, octahe^dron or a.sphere. A co-ordinate system that-is appropriate- 
for the geometry employed simplifies the equations further. The discussion of the method below 
uses a tetrahedron and the formulae are derived in Cartesian co-ordinates. 

Figure 4 shows a tetrahedron with corner points Pi, P2, P3, and P4. We denote the Cartesian co- 
ordinates of point Pi by (x,, yi, zi). Let P.P- denote the vector which starts at Pj and ends at Pj. The 
volume of the tetrahedron in Figure 4. Ktdra, can be calculated from the absolute value of the 



determinant of the matrix M with the vectors PjP, , PjP, , and P^P, as rows: 



^^icira= 7 |det M\. (8) 
o 
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Gauss' theorem states 

-- — — - - -- — - 

\\l{v^u)dy = \[u.c/s (9) 

or alternatively, 

UX(VxuVF = J£uxc/J (10) 

If the d imensions of the tetrahedron configuration in Figure 4 are small compared to the vvaveienmh 
of u, the mean values of the curl and divergence in equations (9) and (10) can be calculated from^ 



"^■"^Tr-Wlc^uw (11) 

tetra 

and 

Vxu--^JJ|^(VxuVK . (12) 

!etra 

Moreover, the surface integrals on the right hand side of equations (9) and (10) reduce to four 
surface integrals respectively, over plane triangular surfaces. 

Since the normal to the surface of the tetrahedron does not vary over each of the triangular sub- 
surfaces, equation (9) becomes 

JJ|(V.u)c/K = J/ udS.n,,, (13) 
+ JJ, ua'5n„3 

where n,^;t is the outward unit normal to the triangular subsurface S,jk with corner points P;, P; and P;; 
(see Figure 4). 

AVhere the side of the tetrahedron is small compared to tlie wavelength of the elastic -vvaveficld, the 
value of u at any point on one of the triangular surfaces of the tetrahedron can be calculated by 
linear interpolation of the wavefield between the three corner points of the triangle. The surface 
integral over 6^- (using linear interpolation between P,, and P^-) then reduces to 



rr A 

tjk 



where A is the area of 5ijk- Since /Inijk = T^P, x 1\P^ /2 we have 
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(,vi + Xi +.V4,>^, +y2 +y4. z\ +Z3 +Z4) • (P3P, x P,P, )] 
Equivalently, the curl of u in P^can be calculated as 

\\[{yxn)dV = ~[ (16) 



+ (xi + -X2+X3,^| +>'2+>'3. ^1 +r2+Z3)x(P,P, X P, P, ) 



+ (X2 + X4 + A-3,>'2 +^4 + + Z4 + Zj) X ( P, P, X P3P2 ) 



+ + A-4+X2,>^| -^>^4+y2.Z, + r^ + Z2)x(P,P( X P, P, ) 
+ (,r, + ,V3+A-4.7l +^3+^4.^1 +Z3 +Z4)X(P;P^ X PJ^ )] 

Finally, the average values of the divergence and curl of u within Kare obtained from equations 
(11) and (12). 

In this way, the divergence and curl of a wavefield can be calculated from a set of closely spaced 
measurements distributed in all spatial directions. A minimum of four measurement points is 
needed. The volume that is bounded by planes connecting these points is generally a tetrahedron 
and where the dimensions of the tetrahedron are small with respect to the wavelength of the 
wavefield, an assumption that the wavefield varies linearly along each side of the enclosed volume 
is possible. The -average divergence and average curl can then be efficiently calculated'using 
Gauss' theorem. 

For instance the average divergence inside the tetrahedron is calculated from the volume integral of 
the divergence of the wavefield over the tetrahedron, divided by its volume. The volume integral 
can in. turn be converted to four surface integrals over each side of the tetrahedron by using Gauss' 
theorem. If the wavefield varies linearly over each side, these surface integrals are simply the 
volume of a prism where the base of the prism is the surface area and the heights of the three 
corners are the magnitude of the wavefield in the respective corners projected onto the normal of 
the surface. 



The average curl of the wavefield is evaluated similarly. The same procedure can be used to 
ca^rculaie"^ve7age~^url~and"d1wrge7i^^^^ - — — — — — — 

Uncertainties in the actual geometry can be eliminated by allowing the exact locations of Pi. P:. P;, 
and P4 to vary subject to either minimising equation (16) over a window(s) with an identified P- 
event and/or minimising equation (15) over S-events. This technique may also serve as an 
alternative means for adjusting for the differences in geophone coupling over the tetrahedron. 

For anisotropic media the decomposition of u into P- and S-wave components as described bv 
equations (2), (3), (4) and (5) is strictly not valid since the (quasi) P- and S-waves in anisotropic 
media both have compressional and shear components. However, this is of secondary importance 
and the theory above will be approximately valid for most practical cases. 

hi heterogeneous media, the heterogeneities act as secondary sources and scatter the wavefield. The 
separation of a wavefield into P-waves and S-waves is not strictly valid in the vicinity of a source. 
However the P/S separation technique described herein is generally applicable for common 
heterogeneous and anisotropic Earth materials, as well as for isotropic materials as discussed above. 

With the apparatus and associated method, P-wave and S-wave components in seismic data can be 
separated. 

In addition attenuation of near-surface scattering is possible as much noise is of S-w-ave character 
and where the primary interest is to record P-wave data, S-waves can be removed by emploving the 
configuration of four one-component geophones as shown in Figure 3. 

Applications of seismic techniques do not only include hydrocarbon exploration but also 
enviromiiental, engineering and groundwater exploration and imaging. For such application's, ic is 
typically of interest to image or detect very shallow structures, as shallow as I- 100m. Near-surface 
generated noise is particularly severe in these applications and a configuration of receivers such as 
in Figure 3 is particularly advantageous for shallow seismics. 

The invention also allows for the attenuation of ground-roll. Since ground-roll can be elliptically 
polarised, it will appear both in the separated P- and S-wave sections, but with lower amplitude in 
the separated P-wave sections. 

P- and S-wave separation is of particular use for marine seismic data recorded at the seafioor. since 
the near-surface at the seafloor is generally is of a more homogeneous character than at land. The 
technique is of use in removing Scholte waves that typically contaminate ocean data. Moreover, 
measurements of the wavefield over a volume also allows for the separation of the wavefield into 
up- and down-going waves. 
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C744/S (57.0307 UK) 

Claims 



1 . Seismic detection apparatus comprising seismic detection means capable of detecting a plurality 
of seismic components over a selected volume. 

2. Seismic detection apparatus according to claim 1, wherein the seismic detection means comprises 
a plurality of receivers spaced apart to enclose the selected volume. 

3. Seismic detection apparatus according to claim 2, wherein the receivers include geophones. 

4. Seismic detection apparatus according to claim 2 or claim 3, wherein the receivers include 
hydrophones. 

5. Seismic detection apparatus according to any of claims 2 to 4, wherein the receivers include 
accelerometers. 

6. Seismic detection apparatus according to any of the preceding claims, wherein the seismic 
detection means is adapted to detect seismic signals along tliree mutually orthogonal axes. 

7. Seismic detection apparatus according to any of claims 2 to 6, wherein the selected volume is 
formed by at least three receivers spaced within a plane and at least one receiver placed external to 
the plane, thereby to enclose the selected volume over which measurement of seismic signals 
occurs. 

8. Seismic detection apparatus according to claim 7, wherein the receivers are spaced to form a 
tetrahedron, with three receivers spaced to define a triangle within the plane and a fourth receiver 
placed external the plane but in line with the centroid of the triangle. 

9. Seismic detection apparatus according to any of claims 2 to 8, further comprising processing 
means which analyses detected seismic components from individual receivers to separate P-wave 
components from S-wave components. 

10. Seismic detection apparatus according to claim 9, wherein the processing means is provided at 
substantially the samejocation as the receivers to allow for.on-site processing. - - 

1 1. Seismic detection apparatus according to claim 9, wherein the processing means is provided 
remote from the receiver to allow off-site processing of detected seismic components. 

12. Seismic detection apparatus according to any of claims 2 to 5 and 7 to 11, wherein the receivers 
are adapted to detect components along one axis and orientation of the receivers is selected so as to 
permit processing of one-component data to identify P-wave components within the plurality of 
seismic components. 

13. Seismic detection apparatus according to any of claims 3 to 12, wherein the spacing of the 
receivers is selected to be smaller than the wavelength of the detected seismic components. 
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14. Seismic detection apparatus according to any of the preceding claims, when attached to a 
wirelin^fot^la^cing^temporar — — - — — — — ^ 

15. Use ofseismic detection apparatus according to any of the preceding claims on surface, 

16. Use of seismic detection apparatus according to any of claims 1 to 13 in a marine environment 
or transition zone, 

17. A method of processing seismic data comprising acquiring seismic data relating to a wavefield 
over a selected volume of acquisition, and measuring the curl and divergence of the wavefield from 
the seismic data, to thereby identify seismic components within the seismic data. 

18. A method of processing seismic data according to claim 1 7, wherein P-wave and S-wave 
components are separately identified. 

19. A method of processing seismic data according to claim 17, wherein up-going and down-going 
wavefield components are identified from the seismic data. 

20. A method of processing seismic data according to claim 17, enabling attenuation of unwanted 
seismic components from seismic data, 

21. A method of processing seismic data according to any of claims 17 to 20, further comprising 
averaging the curl and divergence over the selected volume of acquisition. 

22. A method of processing seismic data according to any of claims 17 to 21, further comprising 
acquiring seismic data over a selected period of time. 

23. A method of processing seismic data according to any of claims 17 to 22, further comprising 
defining the selected volume of acquisition to be small compared to the wavelength of seismic 
waves to be detected. 

24. Apparatus and method substantially as herein described and particularly with reference to and as 
illustrated in the accompanying drawings. 
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Abstract 

Title: Seismic Detection Apparatus and Method 

Seismic detection apparatus comprising seismic detection means capable of detecting a plurality of 
seismic components over a defined tetrahedral volume (36) is provided. The seismic detection 
means comprises four three-component geophones (22, 24, 26, 32), Seismic data acquired by the 
geophones is processed to separate P-vvave components from S-vvave components. The geophones 
are spaced apart by distances smaller than the wavelength of the detected seismic components. The 
apparatus may be used on surface or in a marine environment or transition zone. A method of 
processing seismic data is also provided comprising acquiring seismic data relating to a wavefield 
over a selected volume of acquisition, and measuring the curl and divergence of the wavetleld from 
the seismic data, to thereby identify seismic components within the seismic data. 



[ Figure 2] 
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